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Abstract 

We derive a general theorem relating the energy and momentum with the velocity 
of any solitary wave solution of the generalized KdV equation in iV-dimensions that 
follows from an action principle. Further, we show that our TV-dimensional Lagrangian 
formulation leads to a subclass of the equations discussed recently by Rosenau, Hyman 
and Staley. 
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1 Introduction 



Recently, in an interesting paper, Rosenau, Hyman and Staley fl] have shown the existence of 
spherically symmetric compactons in arbitrary number of dimensions but traveling uniformly 
in x direction. The purpose of this article is to show that a subclass of these compacton 
equations (in arbitrary number of dimensions) can be derived from a Lagrangian. For this 
subclass, we derive a general theorem relating the energy and momentum, with the velocity 
of any solitary wave solution even if the solution is not known in a closed form, and using 
this information comment about the stability of this subclass of compacton solutions. 

Compactons were originally discovered by Rosenau and Hyman (RH) [2 in their investigation 
of a class of generalized KdV equations described by the parameters (m, n). 

ut + {u m ) x + (u n ) xxx = 0. (1) 

They found [21 E] among other things, that these equations have solitary wave solutions of 
the form A[cos(<i^)] 2// ^ m ~ 1 - ) for m — 2,3, where £ = x — ct and —n/2 < d£ < ir/2. However 
these systems of equations were not in general derivable from a Lagrangian. For that reason, 
Cooper, Shepard and Sodano (CSS) jlj were led to consider a related system of equations 
derivable from a Lagrangian. By starting with the first order Lagrangian [3J [S] 

c(i,p) = J ^<p x <Pt + ~ a M p (^*) 2 ) dx = J Ldx ' ( 2 ) 

CSS derived and studied a generalized sequence of KdV equations of the form: 

u t + u l ~ 2 u x + a[2u p u xxx + Apu p ~ l u x u xx + p(p - l)^" 2 ^) 3 ] = , (3) 

where the usual field u(x, t) of the generalized KdV equation is defined by u(x, t) = ip x (x, t). 
The equations of CSS had the same class of solutions as the original RH equations but had 
the advantage of being derived from an action and thus conserved energy. The CSS equations 
lead to similar compacton solutions when the parameters (m, n) of the RH equation have 
the value (I — l,p + 1). In general the RH and CSS equations have different conservation 
laws. 

Recently, Rosenau, Hyman and Staley 1] have shown that a generalization of their original 
set of equations, namely 

C M (m, a + b):u t + (u m ) x + ~[u a (V 2 u% = , (4) 

admitted traveling compacton solutions which were spherically symmetric but traveling uni- 
formly in x direction. Here m, a and b are integers. Introducing the variables: 

s = x — At; R = Js 2 + x\ + ... + x 2 N , (5) 
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one can integrate their equation in a traveling frame to find 



u a 



- Xu +u + m^dR R m u 



0. (6) 



Let us assume that the general solution of this equation is of the form 

u = AU[f3R] , (7) 
and further let the velocity dependence of A and f3 be 

A (x X 7 , 13 X 1 . (8) 

Then one can immediately scale out A provided 

1 m — n 

a= 7 ; r l = 7T { TT > (9) 

m — 1 2{m — 1) 

where n = a + b. It is amusing to note that the velocity dependence of both the amplitude A 
and the width (3 is dimension independent. Further, when m = n the width is independent 
of the velocity of the compacton. (Note that N denotes the number of space dimensions 
which is different from n = a + b). From the RHS Eq. (J3J it immediately follows that one of 
the conserved quantity in the RHS model is the mass M defined by 



M = J ud N x. (10) 

Using Eqs. (jSJ) and Q it immediately follows that the velocity dependence of M is however 
dimension dependent 

M oc \l N ( n - m )+ 2 ]/' 2 ( m - 1 ) . (11) 
Nevertheless, for compactons since m = n, this ^-dependence disappears even for M. 

RHS found two explicit solutions in TV-dimensions: (a) true compacton solution in the case 
a = 1, i.e. Cj^{m = 1 + 6, a + b = 1 + 6) and (b) Cj^(m = 2,n = a + b = 3) when the solution 
is part of a parabola. 

We shall now show that a subclass of RHS Eq. satisfying b = a + 1 can be derived from a 
generalization of the CSS Lagrangian (j2J) with (ip xx ) 2 — > (V(f) x ■ V(b x ). In particular, consider 
the Action 

S = J L(l,p)dxdt = J (^f x <f t + j^ZJj ~ «(^) P (V^ • Vy^ dxdt. (12) 

From this we obtain the equation: 

u t + (d^u 1 - 1 /^ - 1)] + a(d x )[pu v - l Vu ■ Vu] + a(d x )[2u p V 2 u] = . (13) 

Although this is the natural generalization of the CSS equation, at this point to make the 
comparison with RHS more explicit, we will modify the coefficients of the Lagrangian density 
to be: ^ 

L 2 (l,p) = hp x ip t + - a{ip x y{Vip x ■ V Vx ) . (14) 
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The equation of motion that results from this can be written in the form: 



u t + [d^u 1 - 1 ] + ^-^d x y/ 2 V 2 u^ 2+ V] = . (15) 

Thus by identifying 

a = p/2, l = m + l, b = a + l, a = 1/2, (16) 

we obtain a subclass of the RHS equation with b = a + 1. We have not found a way to 
generalize the CSS Lagrangian to obtain the full family of RHS equations with arbitrary b 
which is independent of a. For b = a + 1, the Lagrangian (jUj) was also recently written down 
in the Appendix of j^j. 



2 Conservation Laws 



From Eq. (|15j) we find that 



M 



u a x 



(17) 



is conserved. This conservation law also follows from the invariance of the action under 
if — > ip + <p , where (p is a constant. 

On multiplying Eq. (|15|) by u it is easily shown that another conserved quantity is P = 
| / u 2 d N x. In particular, the corresponding continuity equation is of the form 

Here summation over k is understood with k counting the number of space indices. In this 
case the relevant matter density p and the current density jj~ are given by 



u 



2a 



P+i Q2u 



Jz 



9 2 ' 3y p+l~ dxdy 



-u* 



(p + 1) 



2a 

' p+1' 
2a 

(P + 1) 



-i d 2 u 
dxdz ' 

dx 2 



+ a(p - 1)m p (Vm • Vtz) . (19) 



A third conserved quantity is obviously the Hamiltonian corresponding to the Lagrangian 
density L 2 as given by Eq. (JHJ). It is given by 



H 



d N x 



I 



(20) 



More generically, the Lagrangian density L 2 in Eq. ()14|) is of the form: 



L = L( 



(21) 
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Here we will use the compact notation that \i runs from 0, 1 ... AT with x = t and x\ = x. 
Also M = d(p/dx^, etc. The stationary condition on the action under variations of of the 
generic form 0^0 + 5<p, with <50 vanishing on the boundaries of the integration range, leads 
to the Euler-Lagrange equation for this type of Lagrangian: 

= 0,d^. (22) 

Here we use the Einstein summation convention that repeated indices are summed over. 
To obtain the energy-momentum tensor we consider an infinitesimal arbitrary coordinate 
transformation 

Xp^Xp + e^x). (23) 

Under this transformation 

<t>{x) ^(j)(x + e(x)) =<f)(x) + 5(f)(x), (24) 
which for infinitesimal transformations leads to 

8<f>(x) = ct>^. (25) 



Under this transformation we obtain: 



SS 



d N+1 



x 



-^<9 M (0 A e A ) + -^-dpdvtyxex) 



'\1V 



(26) 



One identifies the energy- momentum tensor T^a from the coefficient of d^ex. That is after 
integration by parts and using the Euler-Lagrange equation one can cast 5S in the form 



5S = J d N+1 x[T^e x } : 



where 



6L 6 L 

T mA = S^L + — A + 2^—0^ - d v 

The A" + 1 conservation laws that arise from the fact that 

dfjT^x = 



5L 



-<t» 



are dP^/dt = where 



P„ 



d N xT 0fl . 



(27) 



(28) 



(29) 



(30) 



The conserved energy H is Pq and there are N conserved linear momenta P\~. In particular 
for the Lagrangian considered in this paper, the only fields entering the Lagrangian are 0o, 
(f> x , and 4> x k- This results in the fact that 



Too = -L+ -ct>At = U 



(31) 
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and 

T ok = \<PAk] k=l ... N. (32) 

Note that ifk, k 7^ 1 is only known through an integral operator, see Eq. (36) below. The 
tensor T ik can be constructed by using 

6L/6<pi = 6 a [<p t /2 + W~7(* - !) - (xp{<PxY- X V<p x ■ V<^ c ] (33) 

and 

SL/5cp iv = -2a5 a 5 uk cp xk (cp x ) p . (34) 

Because the compactons have a preferred x direction, there is also a rotational symmetry 
about the x axis leading to conservation of angular momentum about that axis. 

We shall see below that in general there are only two conserved quantities corresponding to 
the RHS Eq. ©. 



3 General Properties of the Exact Solutions 

Let us assume that the general solution to the field Eq. (JT5j) is of the form 

u = Vx = AZ(PR) , (35) 
where A and j3 are constants, Z is an unspecified function and 

N I 1 / 2 



R 



((* + ?(*)) 2 + £ 



x k 



8=2 



(36) 



Using the fact that (3 must be chosen to minimize the action we shall deduce the general 
behavior of the conserved quantities M, P and H as well as the amplitude A and the width 
(3 as a function of the velocity q. 

To evaluate the action we must first determine (p. A convenient choice for the lower value 
limit of integration (see also leads to: 

<p{x,Xi,t) = f X dx'[AZ{f3R{x' + q{t),Xi)] . (37) 

J-q(t) 

Introducing the rescaled variables y' — (3{x' + q(t)), y% — /3xi, we have 

A r P(x+q(t)) 

¥>(M) = |jf dy Z[tf + Y,vt) 1/2 \. (38) 

Thus we find that 

(ft = q (f x , (39) 
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and the the kinetic term in the action is 

KE = J d N x l -ip xVt = q J d N x~{<p x ) 2 = qP . (40) 

Thus the effective Lagrangian has the canonical form for point mechanics: 

L = Pq-H, (41) 

where 

H = J d N x TC , (42) 

is evaluated for the field ansatz. 

We can now eliminate A in all the conservation equations in favor of the conserved momentum 
P. Defining the constants C\ via 

Q = Jdr Z\r)r N - l tt N , (43) 

where Qn is the solid angle in N dimensions, we find that 

„ A 2 



2(3 N 

so that 



C 2 , (44) 



A = J A^Ar/2 p i/ 2 (45) 



The constant of motion M is given by 



r Pi P 1 / 2 

.1/ J'I X -'-" \ — ( '—- { 



PN/2 ■ 



The key point to note is that (3 can be determined in terms of P using the fact that the 
actual solution of Hamilton's equations satisfies 

This is a consequence of the fact that if we consider a reduced class of trial functions 
parametrized by (3 then the solution that minimizes the energy satifies Eq. (|47jl . When 
the actual solution is in the class of trial functions, then the value of (3 that leads to a solu- 
tion of the actual equations of motion is the value of (3 obtained from the minimization of 
the energy This is due to the fact that the exact equations of motion minimize the Action. 

For the Hamiltonian we find: 

H = -h^PlI 2 ^ 1 ^' 2 + /2 ( p )p( P +2)/ 2/ ^ P+ 4)/2 ; (48) 
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where 



and 



hip) = aD P ( 



2 \ (P+2)/2 



dr r JV " 1 [Z(r)] p (dZ/dryQ N . 



Minimizing if with respect to /3 we obtain: 

/3 = K(N,p, /) 2 /[ Ar ('-P- 2 )- 4 ] p(p-i+2)/[AT(i-p-2)-4] 



where 



Therefore 



where 



H 



K(N,p,l) 



N(l -p- 2) -4 
iV> + 4 



/ 2 (p)(iVp + 4) 

" iV(/-2)A(/) • 

£ (A ^-JV(Z-2)/[JV(Z-p-2)-4] pr 



_ (p + 2)N -l(N -2) 
iV(p + 2 - Z) + 4 
Note that for iV = 1 we get our previous result that 

p + 2 + l 



Using 



6+p-l 



OH 

X = q= dP> 



then we find, in analogy with the CSS equation, that we have 

A = r-, 

where r is now dependent on iV and given by Eq. (J55j) . We also find that 

P OC A [^(P+2-0-H]/[2(I-2)] 

and 



In addition, we find that 



13 oc A a-P-2)/P<i-3)l ; A oc A VK'-a)] . 



M OC A [ Af (P+2-0+2]/[2(«-2)] 



(49) 
(50) 

(51) 
(52) 
(53) 
(54) 

(55) 
(56) 

(57) 

(58) 

(59) 
(60) 



It is worth noting that while the dependence of M and P on A (i.e., q) is dimension-dependent, 
this is not so for A and (3. Note however that in the special case of compactons since / = p+2 
hence, in that case, the dependence of H, P and M on A is also independent of the number 
of dimensions N. On using the correspondence I — m + 1, p — n — 1 = a + b — 1 between 
RHS and our model, it is easy to check that, as expected, the A dependence of A, (3 and M is 
exactly as given by Eqs. © and (fTTf . Note however that without knowing the Lagrangian, 
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one cannot say anything about the velocity dependence of energy H and momentum P which 
is required to study the question about the stability of such solutions. 



Stability of Solutions: Following the analysis in 8J, it is clear that the criterion for linear 
stability is equivalent to the condition 

dP , s 

^>0. (61) 

Using Eq. (p)%j) . it then follows that the compacton solutions of our model (i.e. RHS model 
with b = a + 1) in iV-dimensions are stable provided 2 < I < p + 2 + A/N. However, it is not 
possible to say anything about the other compacton solutions of the RHS model. 



Comparison and Contrast Between General RHS 
Equation and the Lagrangian Subset 



Before finishing this article, we would like to make several remarks about the RHS model £Q 
vis a vis the subset of models defined by our Lagrangian. 



1. While there are N + 2 conserved quantities (excluding angular momentum about the 
x axis) in our model with b = a + 1, (i.e. H, J ud N x and / u 2 d N x), there are in 
general only two known conserved quantities (i.e. / ud N x and J u b ~ a+1 d N x) in the 
RHS model. While the conservation of / ud N x is rather obvious from the RHS field 
Eq. PJ). the conservation of / u b ~ a+1 d N x is not so obvious. However, one can show 
that the continuity equation (|18|) is satisfied in that case with 



P 



u 



b+l-a 



6+1 



a 



Jy 



mu 



(m+6— a) 



+ 



(m + b — a 
(a + b-2 



2b- 1 



u 



26-1 



d 2 u 
dxdy ' 
b 



-u 



Jz 



(26-1) 



2b- 1 

d 2 u 



-u 



26-1 



d 2 u 
dxdz ' 



-u 



(26-2) 



(Vu ■ Vit) + 



(26-1, 
a + b - 1) 



(26 -1) 



dx 2 
u^V 2 u 



(62) 



In the special case when 6 = a+1, then as expected this expression reduces to the one 
given by Eq. (JT§j) provided we make use of the identification as given by Eq. (fT6*|) . 

It may be noted that in case b — a + 1 = then one can show in the RHS model that 
a third conserved quantity is lnw. In particular, in that case the continuity Eq. (|T%|) 
is satisfied with 



P 



lnw 



mu 



Jy = 

(m-l) 



U 



{2b _ 2) dudu 
dx dy 



Jz 



u 



(26-2) 



du du 



1) 



(26-2) 




+ ( 2b ~ 3 ) u (2b-2) 



(Vu • Vu) . (63) 
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Similarly, when b + 1 — a = 1, one can show that a third conserved quantity is u\nu. 
In particular, continuity Eq. (|TH|) is satisfied with 



p = u In u , j y 



dx dy 



,(26-2) 



9m 9m 



^ z dx dz 



-u 



(26-2) 



<9x 



m 
"(26 



+ (b - 1) Inu 



(26-2) 



(Vu • Vu) . 



(64) 



2. In the special case when a = and b = m, one can easily show that the RHS model then 
has two more conserved quantities given by / ucos(y/bx)d N x and / u sm(\/bx)d N x. In 
particular, the corresponding densities are 



p = ucos(Vbu) , j y = bvS b sin(v / &x) , j z = bu^ b sin(Vbx) , ... 



3 x 



-u 



h cos(Vbx) + 7 cos(v / &x)VV + Vbu^" 1 ^ sin(v / 6. 

OT 



(65) 



p = u sin 



j y = fa/ 6 cos(v / 6x) , j 2 = fa/ 6 cos(\/fcc) , ... 



dy 



-u 



b sm(Vbx) + \ sin( Vbx)V 2 u b - VW 6-1 ^ cos(Vbx) . 
b ox 



(66) 



We believe that the existence of these two extra conserved quantities could perhaps 
explain the numerical observation of RHS as to why a = 0, m = n compactons are 
much closer to being elastic compared to those compactons with m = n but a ^ 0. 

3. In the case of the RHS model, the compacton solutions are easily written down in 
iV-dimensions in case a — 1, m — n — b+1. In particular, following RHS pQ, if we 
write down the solution as 



u 



A 



1 - 



F(R*) 



< R < R* 



(67) 



and vanishing elsewhere, then one can show that the solution in iV-dimensions is 

F(R) = R- k J k (VbR) , k = (N — 2)/2 = 0, 1, 2, ... , 
1 d 1 ' 



F(R) 



RdR 



[ci sin(VfaR) + c 2 cos(VbR] , k = (N - 3)/2 = 0, 1, 2, ... (68) 



where J n {x) is the Bessel function of order n [§]• 

4. Similarly, it is easy to generalize the parabolic compacton solution of RHS in the case 
m = 2, n = a + b = 3. In particular, it is easily shown that in general, if m — (n + 1)/2 
then the parabolic compacton solution to Eq. (0J) is given by 



u 



(n-l)/2 = A _ BR 2 ^ < i? < i?, 



A/B, 



(69) 
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where 



.4 



X[N(m-l) + b + l-a] 
6 + 1 - a 



B 



(m - If 



2b[N(m - 1) + b + 1 - a] 



(70) 



5. Apart from these two compacton solutions, it is also possible to obtain other compacton 
solutions to RHS Eq. (jlj). For example, on using the ansatz 



At + a 2 y + a 3 z + 



(71) 



in Eq. (j3J), it is easily shown that the RHS Eq. (j3J) essentially reduces to a one- 
dimensional problem. In particular, on integrating twice, and assuming the two con- 
stants of integration to be zero, one can show that the RHS Eq. fl^J) takes the form 



ds 



2\u 



3—n 



2u 



m+2—n 



b(b + l-a)C b(m + b-a)C 



(72) 



where C = 1 + a\ + a| + .... This is a one-dimensional equation which represents 
the generalization of the old RH equation [2]. Following our earlier work jTUj, one can 
immediately write down the compacton as well as the elliptic and even general hyper- 
elliptic compacton solutions of Eq. (fT2*|) . For example, for m — n, the compacton 
solution is given by 



u = ^[cos^s)] 2 ^"- 1 ) 
and zero elsewhere, where 

, " 26A 
A 



tt/2 < (3s < tt/2 



b+l-a 



V(»-i) 



(n-1) 



2bVC 



On the other hand, for m = 2n — 1, the elliptic compacton solution is given by 

u = A[cn(Ps, k 2 = l/2)} 2/(n - 1) , - K(k 2 = 1/2) < (3s < K(k 2 = 1/2) , 
and zero elsewhere, where 



A 



(36 + a - 1)A 



6 + 1 



l/2(n-l) 



[n-1) 



A 



(6+1 -a)(36 + a- 1) 



1/4 



(73) 



(74) 



(75) 



(76) 



Here cn(x, k) denotes a Jacobi elliptic function with modulus k and K (k) is the com- 
plete elliptic integral of the first kind 0- Similarly, the other solution of ^U] (when 
m = 3n — 2) can be written down as well. 
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5 Conclusion 



In this article we have generalized the Lagrangian of CSS [1] to iV dimensions and have 
obtained a system of KdV equations which are a subset of the recently discovered equations of 
RHS [I] which support iV-dimensional compactons. By exploiting the Lagrangian structure 
we were able to obtain certain general relations among the conserved quantities (derived 
from an energy-momentum tensor) and how these quantities depend on the velocity of a 
traveling wave solution. Using this information we remarked on the question of the stability 
of such solutions, specifically one would expect the a = compactons to be more stable due 
to the existence of two more (N + 2) conserved quantities. 

This work was supported in part by the U.S. Department of Energy. We would like to 
thank J. M. Hyman for sharing his research with us prior to its publication and for fruitful 
discussions. 
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